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Abstract 

Conversion  of  pulses  of  CO2  laser  energy  (18  microsecond  pulses)  to  propellant  kinetic  energy  was  studied 
in  a  Myrabo  Laser  Lightcraft  (MLL)  operating  with  laser  heated  STP  air  and  laser  ablated  delrin  propellants.  The 
MLL  incorporates  an  inverted  parabolic  reflector  that  focuses  laser  energy  into  a  toroidal  volume  where  it  is 
absorbed  by  a  unit  of  propellant  mass  that  is  subsequently  expanded  in  the  geometry  of  the  plug  nozzle  aerospike. 

With  Delrin  propellant,  measurements  of  the  coupling  coefficients  and  the  ablated  mass  as  a  ftmction  of  laser  pulse 
energy  showed  that  the  efficiency  of  conversion  of  laser  energy  to  propellant  kinetic  energy  was  ~  54%.  With  STP 
air,  direct  experimental  measurement  of  efficiency  was  not  possible  because  the  propellant  mass  associated  with 
measured  coupling  coefficients  was  not  known.  Thermodynamics  predicted  that  the  upper  limit  of  the  efficiency  of 
conversion  of  the  internal  energy  of  laser  heated  air  to  jet  kinetic  energy,  a,  is  ~  0.30  for  EQUILIBRIUM  expansion 
to  1  bar  pressure.  For  FROZEN  expansion  a  ~  027.  These  upper  limit  efficiencies  are  nearly  independent  of  the 
initial  specific  energy  from  1  to  1 10  MJ/kg.  With  heating  of  air  at  its  Mach  5  stagnation  density  (5.9  kg/mi>  as 
compared  to  STP  air  density  of  1.18  kg/m°)  these  efficiencies  increase  to  about  0.55  (equilibrium)  and  0.45  (frozen). 

Optimum  blowdown  from  1.18  kg/m3  to  1  bar  occurs  with  expansion  ratios  - 1 .5  to  4  as  internal  energy  increases 
from  1  to  100  MJ/kg.  Optimum  expansion  from  the  higher  density  state  requires  larger  expansion  ratios,  8  to  32. 

Expansion  of  laser  ablated  Delrin  propellant  appears  to  convert  the  absorbed  laser  energy  more  efficiently  to  jet 
kinetic  energy  because  the  effective  density  of  the  ablated  gaseous  Delrin  is  significantly  greater  than  that  of  STP 
air. 

NOMENCLATURE  fin  order  of  usel 

Ef  kinetic  energy  of  vehicle  at  end  of  mission. 
nif  mass  of  vehicle  at  end  of  mission. 

vf  velocity  of  vehicle  at  end  of  mission  in  inertial  frame  of  reference. 

T)  efficiency  of  conversion  of  propellant  kinetic  energy  to  vehicle  kinetic  energy, 
a  efficiency  of  conversion  of  propellant  internal  energy  to  propellant  kinetic  energy, 
ft  efficiency  of  absorption  of  laser  energy  by  propellant 

y  efficiency  of  transmission  of  ground  based  laser  energy  through  atmosphere  to  vehicle. 

El  laser  energy  per  laser  pulse. 

ve  exit  velocity  of  propellant  relative  to  the  rocket,  in  the  rocket  frame  of  reference,  m/s. 
m  mass  of  rocket,  kg. 

F  Force  or  thrust  of  rocket,  N  =  kg  m/s2, 

t  time,  s 

p  propellant  density,  kg/m3. 

<ve>  mass  weighted  average  exit  velocity  of  blowdown  expansion  in  rocket  frame  of  reference,  m/s. 

Ep  kinetic  energy  of  propellant,  J. 

<ve2>  mass  weighted  average  squared  exit  velocity  of  propellant  in  a  blowdown  expansion,  m2/s2. 
mp  mass  of  propellant,  kg. 

I  impulse,  Ns  =  kg  m/s. 

C  coupling  coefficient,  Ns/J  or  s/m. 

O  ratio:  <ve>2/<ve2>.  DISTRIBUTION  STATEMENT  A 

c*  normalized  coupling  coefficient,  c* = c/p.  Approved  for  Public  Release 

V  normalized  exit  velocity,  v = <ve>/$.  Distribution  Unlimited 

(uc-u0)*  normalized  specific  internal  energy,  (uc-u0)*  =  (uc-u°)/  O. 
uc  specific  internal  energy  of  laser  heated  propellant,  J/kg. 

u°  specific  internal  energy  of  propellant  before  laser  heating,  u°,ir  at  STP  =  -  9.0  x  1 04  J/kg. 
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ue  specific  internal  energy  of  propellant  at  the  exit  of  the  rocket  after  isentropic  expansion. 

Vabs  absorption  volume. 

V^*  normalized  absorption  volume,  V2bs*  =  Vabi/p . 

T  temperature,  K. 

P  pressure,  bar. 

h  specific  enthalpy,  J/kg. 

s  specific  entropy,  J/kg  K. 

Mm  average  molecular  weight  of  a  mixture,  g/'mole. 

Cp  specific  heat  capacity  at  constant  pressure. 

va  velocity  of  sound,  m/s. 

X(e-)  mole  fraction  of  electrons, 
e  expansion  ratio 

Ae  area  of  exit  surface 

At  area  of  sonic  surface  or  throat  area 

Subscripts,  Acronyms,  symbols 

i  initial  value  of  property 

f  final  value  of  property 

c  property  in  chamber 

t  property  in  throat 

e  property  in  exit  plane 

p  property  of  propellant 

MLL  Myrabo  Laser  Lightcraft 

SIP  Standard  Temperature  and  Pressure,  298  K,  1 .01 326  bar 
INTRODUCTION 

Laser  propulsion  is  limited  by  laser  power,  so  optimization  of  the  laser  propulsion  mission  may  be  factored 
into  optimization  of  four  energy  conversion  efficiencies,  which,  in  a  first  approximation,  are  independent  of  each 
other.  In  this  idealization  the  kinetic  energy  of  the  propelled  vehicle  at  the  end  of  the  mission  may  be  expressed 
simply: 

(1)  Ef  =  VimfVf2  =  T|aPyEL. 

The  “propulsion  efficiency”,  r[,  is  the  efficiency  with  which  jet  kinetic  energy  is  converted  into  vehicle 
kinetic  energy.  Sutton1  pointed  out,  more  than  50  yean  ago,  that  die  instantaneous  propulsion  efficiency  varies 
during  a  rocket  mission  and  that  it  is  unity  only  when  the  vehicle  velocity  in  the  inertial  frame  is  equal  to  the  jet 
velocity  in  die  rocket  frame.  Unit  propulsion  efficiency  is  achieved  when  die  jet  is  deposited  as  a  stationary  mass 
relative  to  an  observer  in  the  inertia]  frame  of  reference. 

Then,  25-years  ago,  Moeckel2  and  Lo3  independently  and  nearly  simultaneously  published  analyses  of  the 
optimization  of  laser  rocket  propulsion  by  maximizing  q,  and  most  recendy,  Phipps,  Reilly  and  Campbell  (2000, 

2001)4  cited  Moeckel’s  paper  in  their  comprehensive  analysis  of  the  single  stage,  constant  1^  Earth  to  LEO  rocket 
mission.  They  reiterated  the  fundamental  limit  that  Newton’s  second  law  imposes:  for  rocket  missions  that  start  at 
zero  initial  velocity,  the  maximum  q  is  0.648,  which  is  achieved  when  f  =  0.203  and  v/v,  =  1 .595.  For  the  Earth  to 
LEO  mission  the  effective  “delta  v”  (vf)  is  about  10  km/s,  so  the  optimum  single  stage  to  orbit  jet  velocity  is  ~  627 
km/s,  or  specific  impulse  ~  640  s. 

In  this  paper  we  report  a  continuation  of  our  previous  work5  and  analyze  measurements  of  the  overall 
efficiency  of  conversion  of  laser  energy  to  propellant  kinetic  energy,  a|3,  based  on  various  ballistic  pendulum  and 
flight  experiments  with  Myrabo  Laser  Lightcraft,  MLL  [Messitt,  Myrabo,  and  Mead  (2000)6,  Mead,  Squires, 

Beairsto,  and  Thurston  (2000)7].  The  Phipps,  et  al.4  study  defined  an  “ablation  efficiency  “  and  analyzed  the  Earth 
to  LEO  mission  with  unit  ablation  efficiency.  Their  ablation  efficiency  is  equivalent  to  the  product  of  a  and  p. 

Comparison  of  experimental  results  with  thermodynamic  analysis  enables  confining  the  range  of  permissible  p  that  "]  .^.p 
operates  during  the  heating  process.  Wang’ s*  CFD  plasma  model  of  laser  energy  absorption  by  air  have  predicted  j  1 
plasma  temperatures  up  to  30000  K  with  attendant  low  p  values  ~  0.3  to  0.4.  It  has  been  pointed  out  that  p 
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approaches  zero  as  the  plasma  temperature  approaches  -  40,000  K,  where  the  plasma  frequency  approaches  the  laser 
frequency.8, 9 

The  Coupling  Coefficient 

Newton’s  second  law  expresses  the  thrust  that  results  from  expulsion  of  matter  from  a  vehicle  of  mass  m 
at  velocity  vt  as 

(2) 

at 

where  mve  is  the  momentum  of  the  jet  exhaust  in  the  vehicle  frame  of  reference,  [Corliss,  (1 960)]10.  For  the  case 
where  ve  is  constant, 


F  =  -ve- 


Equation  (2)  may  be  used  to  define  an  average  exit  velocity  for  rockets  where  ve  is  not  constant,  such  as 
blowdown  of  a  specified  mass  of  hot  propellant  from  a  fixed  volume,  e.g.,  as  in  laser  rockets  and  pulse  detonation 
rockets: 

t  mf  ff 

JFdt  j  d(mve)  Jd(pvt) 

(4)  <ve>  =  -V  =  JV-  =  AV~- 
Jdm  jdm  jd p 

mi  mi  pj 

so  that  <vt>  is  the  mass  weighted  average  exit  velocity  and  F  =  -  <v„>dm/dt.  Chemical  thermodynamics  may  be 
used  to  rigorously  establish  an  upper  limit  to  <ve>  when  the  propellant  equation  of  state  is  known  and  the  initial  and 
final  states  of  the  propellant  expansion  are  specified. 

The  efficiency  of  conversion  of  laser  energy  to  propellant  kinetic  energy,  ap,  may  be  defined  by  energy 
conservation  in  terms  of  <v?>  for  the  general  case  of  variable  ve 


(5)  Ep  =  ‘/inipCv^  =  ct(5  EL, 

where  the  mass  weighted  average  of  the  square  of  the  propellant  exit  velocity  is 


(6)  <ve2>  =  — 


Jd(pv2) 

a _ 

X  ' 

Pc 


The  impulse,  I  -  jFdt,  imparted  to  a  test  article  by  expansion  of  its  propellant  has  been  accurately  measured 
with  a  ballistic  pendulum  in  the  past5.  Momentum  conservation  requires  equivalence  between  the  measured  impulse 
and  the  propellant  impulse  so  that 

(7)  I  =  mp<ve>. 

Thus,  when  I  and  mp  are  both  measured  an  experimental  <ve>  may  be  obtained 

The  momentum  coupling  coefficient,  also  a  measured  quantity,  is  the  impulse  imparted  to  a  test  article  per 
unit  laser  energy  incident  on  the  propellant. 


(8)  C=~. 
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Using  the  definitions  embodied  in  Equations  (5)  -  (8),  C  may  be  expressed  in  terms  of  a,  B,  <ve>,  and  <ve2>: 


2ap 

■<ve>2  ■ 

_  2ap<t> 

<ve> 

L  <ve2  >  J 

<ve> 

If  ve  is  constant,  <£  =  <\t>2/<\?>  - 1 .  Thermodynamics  may  be  used  to  rigorously  establish  an  upper  limit  to 
d>  and  a  for  any  specified  free-expansion  blowdown  process  when  the  propellant  equation  of  state  is  known5.  The 
$  factor  depends  on  the  distribution  of  exit  velocities,  and  is  mathematically  limited  to  0.5  <  $>  <  1.  It  has  been 
shown51"  that  for  optimum  blowdown  of  laser  heated  air  to  1  bar  pressure  increases  from  0.95  at  low  specific 
energy  (2  MJ/kg)  to  0.98  at  high  specific  energy  (60  MJ/kg).  The  <5  factor  arises  in  Equation  (9)  because  the 
measured  quantity,  the  jet  impulse,  is  proportional  to  mass  weighted  average  velocity  whereas  the  jet  kinetic 
energyis  proportional  to  the  mass  weighted  average  of  the  squared  velocity.  Figure  1  shows  the  relationship 
between  C*  =  C/p,  a,  and  ve*=v,y  $.  For  a  given  a  value,  C*  decreases  as  <ve>*  increases. 

Isentropic  conversion  of  internal  energy  of  propellant  to  kinetic  energy 

Perfect  isentropic  conversion  of  internal  energy  to  propellant  kinetic  energy  occurs  with  no  losses  so  that 

(10)  <ve2>  =  2<uc  -  ue>  =  2a(uc  -  u°),  where 

(11)  a  =  <uc-ue>/(uc  -  u°)  =  <vt2>/2(uc  -  u°). 

These  definitions  generate  a  second  expression  for  C  as  a  function  of  a  that  applies  to  isentropic  energy'  conversion 
from  initial  states  defined  by  (uc  -  u°): 

(12)  C  =  P  [2a<X>/(uc  -  u°)]1/2. 

Figure  1  shows  expansion  isentropes  in  the  C*-a  plane  for  various  initial  state  values  of  (uc  -  u0)*  =  (uc  -  u°)/<I> 
ranging  from  1  MJ/kg  to  1 10  MJ/kg.  At  constant  a,  C*  decreases  as  (uc  -  u0)*  increases. 

Experimental 

Figure  2  shows  a  cross-section  of  the  test  article  (MLL  model  200-3/4)  with  a  ring  of  Delrin  installed  in  the 
shroud.  The  Delrin  shown  weighs  ~  1 0  g,  occupies  a  volume  of  ~  7  cm3,  and  has  a  surface  area  of  ~  25  cm2.  The 
exit  area  of  the  idealized  plug-nozzle11  is  ~  350  cm%  Previous  measurements5  of  Delrin  coupling  coefficients  as  a 
function  of  laser  pulse  energy  (18  jus  pulses)  showed  that  they  rise  to  a  plateau  of  -  350  Ns/MJ  above  EL  -  250  J.  At 
El~  350  J,  the  measured  mass  of  ablated  Delrin  was  ~  40  mg,  which  is  the  mass  of  a  uniform  th in-film  layer  ~  1 1 
micrometers  thick.  Thus,  from  Equation  (7)  ve  ~  3100  m/s  and  from  Equation  (9)  a|3<E>  -  0.54.  If  p  =  1,  C  =  C*, 
and  Figure  1  shows  that  the  initial  specific  internal  energy  that  produces  a  =  0.54/<E>  (with  ®  =  0.98)  in  an  isentropic 
expansion  is  Uc-u°  ~  9  MJ/kg.  At  the  other  extreme,  if  a<t>  =  1,  then  p  =  0.54,  C*  =  1.85  C,  and  Uj-u0  ~  4.8  MJ/kg. 
By  any  analysis,  these  results  show  that  Delrin  is  a  remarkably  efficient  propellant  for  laser  ablation  propulsion. 

For  the  case  of  air,  the  propellant  mass  mp  is  unknown.  Figure  2  may  be  used  to  visualize  a  reasonable 
absorption  volume  for  the  case  where  Delrin  is  absent  and  air  is  the  heated  material.  The  notion  of  an  energy 
absorption  volume,  V^j,  may  be  invoked,  which  contains  a  mass  of  air  propellant  mp  =  pcV^,s  into  which  an  amount 
of  energy  PEl  is  deposited.  In  the  limit  where  the  time  scale  for  energy  absorption  is  much  shorter  than  that  for 
expansion,  the  propellant  density  within  (the  chamber)  remains  constant  during  energy  absorption.  This  enables 
file  initial  specific  internal  energy  of  the  propellant  to  be  defined, 

(13)  Uc  -  u°  =  PEL/pcVlbS) 

where  pc=  1.18  kg/m*  and  u°  =  -  0.09  x  106  J/kg  for  air  at  STP.  Table  1  provides  a  convenient  list  of  values  of  the 
normalized  absorption  volume,  Vabs*  =  V^/p,  derived  from  Equation  ( 1 3)  for  values  of  Uc  -  u°  and  EL  that  lie  within 
the  conceivable  parameter  space  explored  in  experiments.  Table  1  shows  that  an  absorption  volume  for  air,  ~  7  cm*, 
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•would  produce,  with  unit  f5  and  nominal  Et,  values  between  1 00  and  400  J,  heated  air  with  internal  energy  between 
10  and  40  MJ/kg.  If  the  Delrin  surface  shown  in  the  figure,  about  25  cm2*  is  a  suitable  representation  of  the  sonic- 
surface  of  expanding  air,  then,  with  an  idealized  plug-nozzle  exit  area11  of-  350  cm",  the  expansion  ratio  in  this  test 
article  may  be  as  large  as  -  14.  Previously  reported  experiments51  showed  that  air  coupling  coefficients  with  a 
‘loosely”  focused  laser  increased  to  a  plateau  value  of  C<loose  focus)  - 150  Ns/MJ  above  El  - 150  J.  With  a  ■ 
“tightly”  focused  laser  they  increased  to  a  plateau  value  of  C(tight  focus)  -100  Ns/MJ  above  El  -  300  J.  As  shown 
below,  the  maximum  value  of  a$  for  air  in  a  chemical  equilibrium  isentropic  expansion  is  a<J>  -  0.25  and  it  is 
almost  independent  of  initial  specific  internal  energy  and  initial  temperature.  Table  2  summarizes  die  minimum 
values  of  (3,  uc-u°,  C*,  and  temperature  for  various  reasonable  values  of  and  mp.  For  {3  =  1,  mP  =  5.1  mg  (tight 
focus)  or  mp  =  5.6  mg  (loose  focus).  These  are  the  smallest. values  of  mp  consistent  with  the  measured  C  and  El. 

Air  Mollier  Diagram  -  Transformation  of  isentrones  to  the  C*  -  a.  plane 

Figure  3  shows  the  chemical  equilibrium  Mollier  diagram  (u-s  plane)  for  air  up  to  24,000  K.  Figure  3  is 
based  on  the  database  maintained  at  NASA/Glenn  [McBride  and  Gordon  (1995)] 12,  which  is  certified  accurate  up  to 
20,000  K  and  which  is  based  on  extended  9-parameter  fits  to  enthalpy,  heal  capacity,  and  entropy  of  neutral  species 
and  singly  charged  ions.  Above  20,000  K  doubly  charged  ions  begin  to  contribute  but  these  are  not  included  in  the 
database.  This  limitation  leads  to  predictions  of  temperatures  (at  specified  u  and  p)  that  are  too  high  for  plasmas 
above  -  20,000  K. 

Figure  4  shows  a  series  of  isentropes  (vertical  lines)  on  the  Mollier  diagram.  These  are  representations  of 
equilibrium  isentropic  expansions  that  originate  from  initial  states  located  along  the  constant  density  line,  p  =  1.18 
kg/m3,  and  specific  internal  energies  ranging  from  1  to  1 00  MJ/kg.  Table  3  summarizes  other  interesting 
thermodynamic  properties  under  conditions  of  chemical  equilibrium:  T,  P,  h,  s,  Mm,  Cp,  v„  Cp/Cv,  and  X(e).  Table  4 
provides  properties  of  Mach  5  air  at  its  stagnation  density15, 5.9  kg/m3.  Since  the  entropy  of  the  initial  and  final 
states  are  equal,  the  thermodynamic  state  of  the  propellant  in  the  exit  surface  is  uniquely  defined  when  only  one 
additional  property  in  the  exit  surface  is  specified,  such  as  the  exit  pressure  or  the  expansion  ratio,  which  are  also 
indicated  m  Figure  4.  The  expansion  ratio,  e,  is  the  ratio  of  the  area  of  the  exit  surface  to  the  area  of  the  sonic 
surface  or  nozzle  throat,  and  for  isentropic  expansions  this  may  be  represented  in  terms  of  thermodynamic  properties 
in  the  nozzle  throat  and  exit  plane:  z  =  AjAt~  ptVj/ptve . 

,  Figure  5  shows  transformations  of  the  isentropes  in  die  Mollier  plane  to  the  C*-cc  plane  for  equilibrium 
expansions  from  an  initial  density  of  STP  air,  1 . 1 8  kg/mJ.  Lines  of  constant  s  and  pe  are  almost  exactly  coincident 
Lines  of  constant  exit  pressure  run  nearly  parallel  to  lines  of  constant  s  and  pe,  and  all  are  nearly  vertical,  indicating 
that  alpha  is  nearly  independent'  of  ve  and  Uc-u0.  At  five  times  higher  density,  the  constant  exit  pressure  lines  are 
nearly  coincident  with  the  STP  constant  pressure  lines  but  their  exit  pressures  are  about  five  times  higher.  Thus,  in 
the  C*-a  plane,  the  0.2  atm  exit  pressure  line  at  pc  =  1 .1 8  kg/tn3  is  almost  coincident  with  a  1  atm  exit  pressure  line 
at  pc = 5.90  kg/m3.  The  effect  of  higher  pe  is  to  increase  a  from  -  0.4  to  0.6. 

DISCUSSION 

Coupling  coefficients  measured  with  the  Figure  2  test  model  were  reported  in  our  previous  papers5.  With 
increasing  laser  energy  they  rise  to  plateaus  above  about  300  J.  At  El  -  350  J,  C(Delrin)  -  350  Ns/MJ  and  the 
ablated/vaporized  mass  was  mp  -  40  mg.  This  means  feat  <ve>  -  3 1 00  m/s  by  Equations  (7)  and  (8),  and  a[3<J>  - 
0.54  by  Equation  (9).  Thus,  since  (3<$  <  1,  a  >  0.54,  which  is  remarkably  high.  Air  and  Delrin  will  show  very 
.similar  expansion  behavior.  As  shown  previously51*,  fee  dependence  of  a  on  the  density  of  fee  heated  air  is  quite 
strong.  At  [uc-u°]*  =  10  MJ/kg  and  expansion  to  1  bar,  fee  instantaneous  a  increases  from  about  0.43  to  about  0.60 
when  fee  density  increases  from  fee  STP  value  (1.18  kg/m3)  to  its  Mach  5  stapation  value  (5.9  kg/m3).  These 
instantaneous  a  values  decrease  to  about  032  and  030  respectively  for  fee  free-expansion  blowdown  process511.  An 
a  value  around  0.5  is  reasonable  when  the  density  of  fee  ablated  and  vaporized  Delrin  is  as  high  as  -  6  kg/m3  and 
the  blowdown  expansion  is  near  perfect  with  e  - 16.  Most  importantly,  it  appears  that  most  of  the  inefficiency  in 
the  composite  a{33>  efficiency  is  carried  by  a<I>  and  that  3  is  very  close  to  unity. 

Coupling  coefficients  for  air  were  found  to  depend  strongly  on  the  quality  of  the  laser  beam,  as  between  a 
tightly  focused  beam  that  produced  a  lower  C  - 1 00  Ns/MJ  than  a  loosely  focused  beam,  which  produced  C  - 1 50 
Ns/MJ.  This  may  be  due  to  fee  tight  beam  heating  a  smaller  mass  of  air  to  a  higher  energy  than  fee  more  diffuse 
loosely  focused  beam.  Although  the  exit  velocity  would  be  higher  in  the  tight  beam  case,  the  total  impulse  may  be 


lower  because  the  heated  mass  is  lower.  Figure  2  shows  the  geometry  and  size  relationship  of  a  7  cm5  absorption 
volume  inside  the  shroud,  which  contains  -  8  mg  of  air.  With-C(air,  loose  focus)  =  150  Ns/MJ  at  El  -  300  J,  and  tnp 
=  8  mg  we  may  deduce  <vt>  -  5600  m/s,  and  ap$  -  0.42.  If  the  absorption  volume  is  double,  then  <Ve>  and  ap$ 
are  halved.  Figure  5  shows  that  STP  air  heated  to  10  MJ/kg  for  example  would  blowdown  to  1  bar  with  a  =  0.32 
(equilibrium  expansion)  or  a  =  0.27  (frozen  expansion). 

If  we  accept  that  a  reasonable  upper  limit  operational  alpha  is  ~  0.30  in  our  experiments  tx  <  0.j,  then  the 
measured  C  ~  150  Ns/MJ  and  Equations  (7-9)  with  BO  =  1  require  <ve>  <  4000  m/s,  and  mp  >  1 1  mg.  Now  if  BO  is 
~  0.3  as  has  been  suggested  by  CFD  modeling,8  then  the  upper  limit  of  <ve>  decreases  to  1200  m/s  and  the  lower 
limit  of  mp  increases  to  36  mg.  It  would  seem  apparent  that  the  value  of  B  is  somewhat  larger  than  0.3  because  both 
the  upper  limit  <Ve>  and  lower  limit  mp  with  B  =  0.3  are  not  reasonable  for  the  geometry  shown  in  Figure  2. 

CONCLUSIONS 

Experimental  studies  of  the  propulsion  of  a  200-3/4  model  Myrabo  Laser  Lightcraft  healed  by 
10.6  (i  radiation  from  a  C02  laser  showed  that  the  efficiency  of  conversion  of  laser  energy  to  propellant  kinetic 
energy  was  54%  for  Dehin  propellant 

Thermodynamic  analysis  of  isentropic  expansion  of  a  unit  mass  of  air  after  laser  heating  at  constant  volume 
was  examined  under  conditions  where  either  chemical  equilibrium  or  frozen  composition  was  maintained.  The 
upper  limit  for  the  efficiency  of  conversion  of  internal  to  kinetic  energy  in  optimum  blowdown  to  1  bar  pressure  is  a 
-  025,  which  is  nearly  independent  of  the  initial  energy.  For  STP  air,  blowdown  to  1  bar  is  achieved  with 
expansion  ratios  from  s  ~  4  at  high  energy  to  £  ~  8  at  low  energy.  With  this  small  effective  s  ~  4,  equilibrium 
expansion  was  only  slightly  more  efficient  than  frozen  expansion.  Heating  of  propellant  to  higher  energy  states 
resulted  in  only  slightly  lower  a  but  much  higher  exit  velocity.  The  thermodynamic  limitations  were  illustrated  by 
process  representations  of  blowdown  in  the  Mollier  plane. 
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Table  1 .  Normalized  absorption  volume  for  air  at  1. 1 8  kg/m"  as  a 
function  of  internal  energy  and  laser  energy. 


u  _ Vabs/3,  normalized  absorption  volume,  cmJ _ 

MJ/kg  El=50J  El=100J  El=1  50  J  EL=200J  EL=300J  EL^400J 


1 

42.3 

84.7 

127.1 

169.4 

254.2 

338.9 

2 

21.1 

42.3 

63.5 

84.7 

127.1 

169.4 

3 

14.1 

28.2 

42.3 

56.5 

84.7 

112.9 

4 

10.5 

21.1 

31.7 

42.3 

63.5 

84.7 

5 

8.47 

16.9 

25.4 

33.9 

50.8 

67.8 

6 

7.06 

14.1 

21.1 

28.2 

42.3 

56.5 

7 

6.05 

12.1 

18.1 

24.2 

36.3 

48.4 

8 

5.30 

10.5 

15.8 

22.1 

31.7 

42.3 

9 

4.71 

9.42 

14.1 

18.8 

28.2 

37.6 

10 

4.24 

8.47 

12.7 

16.9 

25.4 

33.9 

15 

2.82 

5.65 

8.47 

11.3 

16.9 

22.6 

20 

2.12 

4.24 

6.36 

8.47 

12.7 

16.9 

30 

1.41 

2.82 

4.24 

5.65 

8.47 

11.3 

40 

1.06 

2.12 

3.18 

4.24 

6.36 

8.47 

50 

0.85 

1.69 

2.54 

3.39 

5.08 

6.78 

60 

0.71 

1.41 

2.12 

2.82 

4.24 

5.65 

70 

0.61 

1.21 

1.82 

2.42 

3.63 

4.84 

80 

0.53 

1.06 

1.59 

2.12 

3.18 

4.24 

90 

0.47 

0.94 

1.41 

1.88 

2.82 

3.77 

100 

0.42 

0.85 

1.27 

1.69 

2.54 

3.39 

110 

0.39 

0.77 

1.16 

1.54 

2.31 

3.08 

Table  2.  Measured  and  calculated  quantities  for  expansion  of 
laser  heated  STP  air. 


Quantity 

tight  focus 

loose  focus 

Measured  Quantities 

C(Ns/MJ) 

100 

150 

El(J) 

300 

150 

I  (Ns) 

0.030 

0.022 

Calculated  Quantities  with  (a<5)-... 

=  0.25 

V*ts  (cm3) 

5.08 

5.59 

mp  =pcVabs  (mg) 

6.00 

6.60 

<ve>  =  I/mp  (km/s) 

5.00 

3.33 

Pmir.=  C<V£>/2(a$)m3x 

1.00 

1.00 

C*w  C/p  (Ns/MJ) 

100 

150 

(vu°)min=  pEL/mp  (MJ/kg) 

50 

22.7 

Tmill  see  Table  3  (K) 

14400 

8700 

Vafe  (cm3) 

6.78 

6.78 

mp  =  Pc  Vabs  (mg) 

8.00 

8.00 

<ve>  =  I/mp  (km/s) 

3.75 

2.81 

Pmin  =  C<Vj>/2(a4>)m:K 

0.76 

0.84 

CW-C/p  (Ns/MJ) 

131 

178 

(u=-u°)mm=  PEi7mp  (MJ/kg) 

28.5 

15.7 

Trom  see  Table  3  (K) 

9500 

7500 

Vji,s  (cm3) 

13.6 

13.6 

mp  =  Pc  Vate  (mg) 

16 

16 

✓ — \ 

jf 

II 

£ 

1.8 

1.4 

Pnrni-  C<Ve>/2(ct4>)max 

0.38 

0.42 

C*n«oc=C/p  (Ns/MJ) 

262 

356 

(vn°)min=  PEi/mp  (MJ/kg) 

14.3 

7.9 

T™,  see  Table  3  (K) 

7000 

5400 

Table  3.  Thermodynamic  properties  of  equilibrium  air,  p  =  1.18  kg/rrri. 


u 

T 

P 

h 

s 

Mm 

X(e0 

Va 

Cp/Cv 

MJ/kg 

103K 

bar 

MJ/kg 

KJ/kgK 

KJ/kgK 

kg/kmol 

km/s 

-0.9 

0.298 

1.000 

-.004 

6.864 

1.005 

28.965 

0 

0.346 

1.40 

1 

1.6 

5.4 

1.5 

8.2 

125 

29.0 

4E-10 

0.77 

1.30 

2 

2.5 

8.6 

2.7 

8.7 

1.51 

28.9 

3.E-09 

0.95 

124 

3 

3.2 

11.1 

3.9 

9.0 

2.16 

28.6 

3.E-08 

1.06 

120 

4 

3.7 

13.1 

5.1 

9.3 

2.83 

27.8 

3.E-07 

1.15 

1.19 

5 

4.1 

15.0 

6.3 

9.6 

3.15 

26.9 

2.E-06 

1.23 

1.19 

6 

4.5 

16.9 

7.4 

9.8 

3.04 

26.1 

5.E-06 

1.32 

121 

7 

4.9 

19.1 

8.6 

10.0 

2.69 

25.3 

2.E-05 

1.41 

1.23 

8 

5.4 

21.5 

9.8 

10.2 

2.56 

24.7 

4.E-05 

1.50 

123 

9 

5.9 

23.9 

11.0 

10.4 

2.86 

24.2 

8.E-05 

1.57 

1.21 

10 

6.3 

26.0 

12.2 

10.6 

3.43 

23.8 

l.E-04 

1.62 

1.19 

15 

7.5 

34.1 

17.9 

11.3 

6.70 

21.7 

5.E-04 

1.84 

1.17 

20 

8.3 

41.3 

23.5 

11.9 

8.93 

19.8 

9.E-04 

2.02 

1.17 

30 

9.7 

56.2 

34.8 

13.0 

9.09 

16.9 

3.E-03 

2.38 

1.19 

40 

11.5 

75.4 

46.4 

14.0 

5.13 

15.0 

l.E-02 

2.81 

124 

50 

14.4 

101.0 

58.5 

14.8 

4.81 

14.0 

4.E-02 

3.26 

125 

60 

16.6 

124.0 

70.5 

15.4 

6.62 

132 

l.E-01 

3.60 

124 

70 

18.4 

145.0 

82.3 

16.0 

825 

12.4 

l.E-01 

3.91 

124 

80 

19.9 

167.0 

94.1 

16.5 

9.51 

11.7 

2.E-01 

4.20 

124 

90 

21.3 

189.0 

106.0 

17.0 

10.40 

11.1 

2.E-01 

4.48 

1.25 

100 

22.6 

211.0 

118.0 

17.4 

10.90 

10.5 

3.E-01 

4.76 

1.26 

no 

23.9 

235.0 

130.0 

17.9 

11.10 

10.0 

3.E-01 

5.03 

1.27 

Table  4.  Thermodynamic  properties  of  Mach  5  air  at  stagnation  density,  p 

=  5.90  kg/m3 

.  ' 

u 

T 

P 

h 

s 

M 

X(e) 

Va 

Cp/Cy 

MJ/kg 

103K 

bar 

MJ/kg 

KJ/kgK 

KJ/kgK 

kg/kmol 

km/s 

0.102 

0.560 

9.492 

0.263 

6.864 

1.042 

28.965 

0 

0.471 

1.38 

1 

1.6 

27.1 

1.5 

7.7 

125 

28.97 

4e-13 

0.77 

1.30 

2 

2.6 

43.2 

2.7 

82 

3.45 

28.95 

6.E-1 1 

0.96 

1.25 

3 

3.3 

56.5 

4.0 

8.6 

1.85 

28.73 

2.E-08 

1.08 

1.21 

4 

3.9 

67.7 

5.1 

8.9 

2.33 

28.19 

3.E-07 

1.17 

120 

5 

4.4 

7S2 

6.3 

9.1 

2.65 

27.46 

2.E-06 

126 

1.20 

6 

4.8 

88.9 

7.5 

93 

2.71 

26.69 

6.E-06 

135 

122 

7 

5.3 

100.3 

8.7 

9.5 

2.61 

25.96 

2.E-05 

1.45 

123 

8 

5.8 

112.4 

9.9 

9.7 

2.55 

25.32 

4.E-05 

1.53 

123 

9 

6.3 

124.5 

11.1 

9.9 

2.69 

24.79 

8.E-05 

1.61 

122 

10 

6.7 

135.8 

12.3 

10.0 

3.04 

24.32 

l.E-04 

1.67 

1.21 

15 

8.2 

182.0 

18.1 

30.7 

5.49 

22.19 

6.E-04 

1.91 

1.18 

20 

92 

222.3 

23.8 

31.2 

7.36 

20.32 

l.E-03 

2.11 

1.18 

30 

10.8 

304.9 

35.2 

12.2 

8.05 

17.41 

3  £-03 

2.49 

120 

40 

12.7 

404.9  ' 

46.9 

13.1 

5.52 

15.45 

l.E-02 

2.92 

124 

50 

15.6 

534.8 

59.1 

13.8 

4.28 

14.33 

3.E-02 

3.39 

127 

60 

18.4 

667.9 

71.3 

14.4 

5.20 

13.54 

8.E-02 

3.78 

1.26 

70 

20.8 

794.6 

83.5 

14.9 

6.32 

12.81 

l.E-01 

4.13 

1.27 

80 

22.8 

919.9 

95.6 

15.4 

7.26 

12.14 

2.E-01 

4.45 

1.27 

90 

24.6 

1046.6 

107.7 

15.8 

7.99 

11.52 

2.E-01 

4.76 

128 

Figure  1.  Defined  relationships  between  six  variables  of  interest:  C*,  a,  j3,  &,<vp>*.  and  fu.-u"J*.  The 
plots  show  C*  =  C/p  as  a  function  of  a,  with  lines  of  constant  ve*  =  <vp>/&  and  constant  fuc-ue]*  =  [uc- 
u"]/0.  The  plots  may  alternatively  be  interpreted  as  a  C  vs  a  plots  with  lines  of  constant  vc*  =  <v^>/f<t>  and 
constant [uc-u0]*  =  \uc-u]/ $</>.  .  . 


Figure  2.  Cross-sectional  view  of  Myrabo  Laser  Lighicra.fi,  Model  200-3/4.  The  maximum  diameter  of  the 
test  article  at  the  shroud  is  ~  10  cm.  The  indicated  ring  of  Delrin  weighs  ~  10  g  and  has  a  volume  of  ~  7 
cm 3  and  a  surface  area  ~  25  cm2.  The  idealized  maximum  plug  nozzle  exit  area  is  -  350  err, i2. 
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Figure  3.  Mollier  diagram  for  air  including  singly  ionized  species.  Molecular  weights  are  indicated  at  intersections  of 
isobars  and  isotherms.  The  lower  diagram  shows  a  heavy  constant  density  line,  p  =  1.18  kg/m3  above  a  heavy  constant 
pressure  line,  P  =  1  atm.  The  maximum  energy  initial  states  of  laser  heated  STP  air  lie  on  the  constant  density  line  and 
the  optimally  expanded  states  lie  vertically  below  on  the  constant  pressure  line. 
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Figure  5.  Isentropic  expansions  of  laser  heated  STP  cur  (1.1S  kg/n?) from  initial  states  of  specific 
internal  energy  ranging  from  1  to  110  MJ/kg.  The  circles  and  nearby  crosses  represent  the 
blowdown  quantities  obtained from  initial  [uc-u°J*  states  of  2E3,  and  1 E4  J/kg  for  the  frozen 
expansion  and  2E3, 6E3, 1E4,  and  4E4  kJ/kgfor  the  equilibrium  expansion. 


